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This paper presents a euekoo seareh (CS) algorithm for determining the 
global maximum power point (GMPP) traeking of photovoltaie (PV) under 
partial shading eonditions (PSC). The eonventional methods are failed to 
traek the GMPP under PSC, whieh deerease the reliability of the power 
system and inerease the system losses. The performanee of the CS algorithm 
is eompared with perturb and observe (P&O) algorithm for different eases of 
operations of PV panels under PSC. The CS algorithm is used in this work to 
eontrol direetly the duty eyele of the DC-DC eonverter without proportional 
integral derivative (PID) eontroller. The proposed CS model ean traek the 
GMPP very aeeurate with high effieieney in less time under different 
eonditions as well as in PSC. 
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1. INTRODUCTION 

It is known that the world's electricity consumption will increase annually it is expected that further 
growth will be due to the increase in the number of population and the increasing demands of the modem 
way of lifestyle [1,2]. Therefore, there is a need to develop renewable energy sources to ensure sustainability 
high power supplies to the consumer, as well as to reduce local and global room pollution. The main 
problems in the practical use of Pa V power generation system (PGS) are the relatively low efficiency of 
conversion of primary energy (from 9 to 17%) and the pronounced dependence of the solar cell's energy 
characteristics on external climatic conditions [3-5]. Solar panels consist of several series-parallel connected 
solar modules and are designed for direct conversion of solar energy into electrical energy with the necessary 
values of constant voltage and current. 

The power characteristics of solar cells have a nonlinear convex shape with an extremum, so to 
obtain the maximum available power it is advisable to operate them at the maximum power point (MPP). On 
the other hand, the generated power of solar panels depends on the operating conditions: the intensity of solar 
radiation, ambient temperature [6]. Changing in one or all of these conditions will lead to change in the 
position of MPP. So the search and tracking of MPP in real operating conditions are one of the ways to 
improve the operational efficiency of the solar panel. With no Shading, many conventional methods of 
having been developed for MPPT, the most common methods among them are a constant voltage MPPT, 
perturb and observe (P&O) and its numerous modifications, and incremental inductance method [7-11]. 
Under partial shading, the P-V characteristic has many local peaks and one global peak [12-14]. The 
conventional methods are failed to track the GMPP. Recently, the researchers have worked how to track the 
GMPP under PSC using soft computing techniques. 
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Among them are the fuzzy logie eontrol (FLC) [15], artifieial neural networks (ANN) [16], Genetie 
algorithm (GA) [17], partiele swarm optimizations (PSO) [18] and ant eolony optimizations [19]. They test 
the performanee of their algorithms under dynamie ehanges of the solar irradianee of the PV panel. Although 
the very high eost of implementation the ANN and fuzzy logie, they eapable of traeking the GMPP with 
lower oseillation in power as eompared to P&O algorithm. The major differenee between the proposed CS 
and other teehniques is that the proposed algorithm is used to eontrol the duty eyele direetly without using the 
without Proportional Integral Derivative (PID) eontroller [20]. 

In this paper, the CS algorithm is introdueed to traek the GMPP of the PV under fast variation of 
solar irradianee as well as partial shading eondition. The proposed algorithm is eompared with the P&O 
algorithm in different test eases and the effieieney of both of them are measured in eaeh test ease. Seetion 
two deseribes the modeling of single diode PV eell. The effeet of PSC on PGS is explained in seetion three. 
Seetion four deseribes the CS and P&O and its applieations in PSC. Simulation results of the developed 
model are explained in the sixth seetion. 


2. SINGLE DIODE PV MODEL 

A single diode equivalent eireuit is the most eommon mathematieal model of PV arrays that 
available in the literature. Figure 1 illustrates the equivalent eireuit of the PV eell that used in the 
proposed model. 
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Figure 1. The equivalent eireuit of a solar eell 


The output eurrent of the PV eell is ealeulated by the following equations [21]: 

^=4-4S 
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Where, V, I - is the value of output voltage and eurrent of the PV eell; 4- reverse saturation eurrent of the 
diode; Rs h Rp series and parallel resistanee of the solar eell; n is the ideality faetor; T- is the absolute 
temperature of the PV eell; k=1.38T0'^^ j/k Boltzmann eonstant. 

A PV module eonsists of a large number of identieal solar eells eonneeted in series and in parallel, 
whieh allows inereasing its voltage and output power. The main eleetrieal parameters of the PV module that 
used in this work are summarized in Table 1. 


Table 1. Speeifieations of PV Module 


Parameter 

Value 

Maximum power rating Pmax 

249 w 

Rated voltage Vmpp 

30V 

Open eireuit voltage Voc 

8.3A 

Short eireuit eurrent Isc 

8.83A 

Rated eurrent Impp 

36.8V 

Temperature eoeffieient of Tc 

(0.065±0.015)%/°C 

Temperature eoeffieient of power 

-(0.5±0.05) % / °C 
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3. PV MODULE UNDER PSC 

In real operations, the PV modules are often subject with PSC, which resulted from clouds; trees; 
buildings; towers. As an example to show the effect of the shading on the PV module, four PV panels are 
connected in series. This configuration is named “4sIp” connection, as shown in Figure 2. To simulate that 
model in Matlab, PV output current is supplied to DC voltage controlled source as shown in Figure 3. The P- 
V and V-I for these cases are shown in Figure 4. 
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Figure 2. Normal Irradiance and PSC of PV panels Figure 3. Interfacing mathematical PV module 

model to physical ports 




Figure 4. The P-V and V-I characteristics of PV panels under PSC 


The current generated by PV under PSC is not constant. In case of normal irradiance of the PV 
panels, a constant current of approximately 8.3 A is generated at a functional operating voltage of 0 V to 120 
V as shown in Figure 4. However, when the PV is under the PSC, the generating current cannot maintain a 
constant value. For example, in PSC, PV generates constant current 8.3 A from 0 V to 25 V, then decreases 
to 8 A, 4.1 A and 1.6 A at 55 V, 100 and 120 V, as shown in Figure 4. The PSC case demonstrate four peaks 
because it receives four different levels of irradiation. 


4. CUCKOO SEARCH (CS) 

Cuckoo Search (CS) was first developed by [Yang &Deb] [22]. These algorithms usually work 
based on a random search, which simulates the reproduction strategy of cuckoo birds. It is observed that 
several species of cuckoos perform brood parasitism, i.e. by laying their eggs in host birds’ nests. Usually, 
three types of brood parasitism are seen (1) intraspecific, (2) cooperative and (3) nest takeover. Some host 
birds engage in direct conflict with the intruding cuckoos. If a host bird discovers the cuckoo eggs it will 
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throw these alien eggs away or abandons its nest and builds a new nest elsewhere. Some euekoo speeies have 
evolved in sueh a way that female parasitie euekoos mimie the shape and eolor of the eggs of a few ehosen 
host speeies. To inereases its reproduetion probability. 

4.1. Levy flight 

Searehing for a suitable host bird’s nest is an important part of euekoo’s reproduetion strategy. 
Normally, the seareh for the nest is similar to the seareh for food, whieh takes plaee in a random or in a 
quasi-random form. In general, while searehing for food, animals ehoose direetions or trajeetories that ean be 
modeled by eertain mathematieal funetions. One of the most eommon models is the Levy flight. A reeent 
study by Reynolds and Frye shows that a Levy flight ean be thought of as a random walk where the step size 
has a Levy probability distribution. In CS, nest searehing is eharaeterized by Levy flight. Mathematieally, a 
Levy flight is a random walk where step sizes are extraeted from Levy distribution aeeording to a power law 
as shown below [23]: 

y = i~' ( 3 ) 

Where / is the flight length and X is the varianee. Sinee 1 < A- < 3, y has an infinite varianee. Figure 5 shows 
an example of Levy flight in a two-dimensional plane. 



Figure 5. Example of Levy flight in a two-dimensional plane 


4.2. Cuckoo search algorithm 

Three idealized rules for CS based on euekoo’s brood parasitie behavior: (1) Eaeh euekoo lays one 
egg at a time and plaees it in a randomly ehosen nest (2) The best nest with the highest quality of eggs will be 
earried over to the next generation (3) The number of available nests is fixed and the number of euekoo eggs 
diseovered by the host bird maintains a probability Pa, where 0 < Pa < 1. When generating a new solution for 
a euekoo, a Eevy flight is performed as in the following expression [24] [25]: 

= x/ +a® Levy(A) (4) 

Where x/ are samples/eggs, i is the sample number, t is the number of iteration and a is the step size, whieh 

is related to the seales of the problem of interests. In most eases, we ean use a=l. The produet 0 means entry 
wise multiplieations. The value of a is ealeulated by the following equation: 

a = (5) 

The value of Eevy (>1) is found from Eevy distribution given in Eq. 7 

Zevy (A.) ^u = r^ (6) 


Int J Pow Elee & Dri Syst, Vol. 10, No. 2, June 2019 : 1081 - 1089 





Int J Pow Elec & Dri Syst 


ISSN: 2088-8694 


□ 1085 


5. MPPT USING CS 

All the voltages are provided to the PV panels and the output power are calculated by the proposed 
model. The fitness function in CS is considered to be the output duty cycle that will fed to the De-DC boost 
converter at the GMPP of the PV panels, the dbest is the best one among the fitness. The duty cycle is changed 
in CS by the following equation: 

= dbest + « ® Levy{X) (7) 


Where, di=di,d2...n 

A simplified scheme of the Levy distribution can be explained by: 

■^ = «o i^bes, -^i)® Levy(A) ^kx-d.) 

(|vh 

Where, n=4, y^=1.5 and A:=0.8 is the Levy multiplying coefficient (chosen by the designer), while u 
and V are determined from the normal distribution curves. 

The flow chart of CS that illustrates how to use the algorithm to get the MPP of solar panel is shown 
in Figure 6. 



Figure 6. The flowchart of the CS algorithm 
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6. SIMULATIONS RESULTS 

The proposed model is earried out in the Matlab/Simulink program software as shown in Figure 7. 
The DC-DC boost eonverter is seleeted in the proposed work to mateh the MPP eontroller with a resistive 
load. The seleeted parameters of the eonverter are; 

L=l.l mH; ei=0.4*10 '^ F, e2=0.4mF and switehing frequeney f=25 kHz. The eriteria for stopping the 
iterative proeess of ealeulations are the maximum number of iterations and the relative error in ealeulating 
the global extremum of the objeetive funetion: 

-Pbes, .100 <0 3 (10) 

^best 


Figure 8 and 9 show the performanees of the CS and P&O for determination the GMPP of the PV 
under normal irradianee of the solar panels and in ease of PSC that was shown in Figure 2. The following 
equation is used for evaluating the effieieney of the two algorithms in all eases 

/ 7 =^- 100 , % ( 11 ) 

max 

Where Pq is the value of the output power of the PV panel, that traeked by the algorithm; Pmax - the 
value of the maximum available power. 
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Figure 7. MATLAB/Simulink of the proposed model 


As illustrated in Figure 8 the CS algorithm is able to traek the GMPP of the solar panel after only 
four iterations the optimal duty eyele was (dbest=0.47) and the traeking time was 0.32 see. Also in ease of PSC 
the proposed algorithm traeked the global peak, while the P&O algorithm trapped in loeal peak whieh 
inerease the system losses and reduee the reliability of the PV plant. 


(a) 


(b) 


Figure 8. The performanee of the two algorithms in the normal ease (a) CS (b) P&O 
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(a) (b) 

Figure 9. The performance of the two algorithms in the PSC (a) CS (b) P&O 


6.1. Fast variation of the solar irradiance under constant temperature 

Case 3 are about the fast variation of the solar irradiance of the PV panel under PSC. In case 3, as 
illustrated in Figure 10. From time t=0 to t=l sec all PV panels received the same amount of the solar 
irradiance G=1000w/m2, after t=l sec the PSC occurs where, the solar irradiance of the PV panels were 
Gl=1000 w/m2, G2=500 w/m2, G3=1000 w/m2 and G4=100 w/m. On the other hand, in case 3 as shown in 
Figure 10 after shifting from normal operations to PSC, the CS only can handle the GMPP in this case with 
highly tracking efficiency (/7=99.7%), while the conventional P&O trapped at the local peak. Also, can be 
seen, the steady power performance of the CS are very well. Table 2 shows illustrates the detailed 
comparison between the CS and P&O for all test cases. 



fcA- - 
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(a) (b) 

Figure 10. (a) The performance of CS and P&O under fast variation of the PV under PSC. (b) duty cycle of 

the CS algorithm 


Table 2. Results of the three cases of studies 




CS 



P&O 


case studies 

Pmax, 

Watt 

Po, 

Watt 


Pmax, 

Watt 

Po, 

Watt 


Normal Irradianee 

996 

995 

99.8 

996 

976 

97.9 

PSC 

452 

451 

99.7 

452 

331 

73.9 

Dynamic operation with PSC 

484 

483 

99.7 

484 

334 

69 


7. CONCLUSIONS 

In this study, the performance of CS algorism investigated to track the GMPP of the solar panel 
under PSC and compared with the conventional P&O algorithm. Three different test cases of studies of PV 
modules are used under PSC to evaluate the performance of the proposed CS algorithm. For all test cases, the 
tracking efficiency of the CS is higher than 99 % within 300-500ms. The developed CS algorithm can handle 
the PSC very efficiently under different shading conditions, while P&O algorithm incapable of handling the 
PSC of the solar panel. The results also indicate that the tracked power by the CS has a very low fluctuations 
of the steady power compared to P&O algorithm. The CS has proven accuracy, robustness, and effectiveness 
of efficient energy utilization for the standalone PV system. 
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